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New insights on chemical evolution of galaxies: XMM-Newton 
observations of M82 
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We report on the ongoing data analysis of a very deep (~ 100 ks) XMM-Newton observation of the starburst 
galaxy M82. We show some details of data analysis and a few results from spatially-resolved spectroscopy with 
the EPIC cameras. Since M82 is a bright object with a complex spectrum, the data reduction of such a deep 
observation has posed many challenges both about the involved astrophysical processes and the available data 
analysis techniques. Vertical (with respect to the galaxy plane) abundance gradients are discovered. The hints 
for an under-abundance of Oxygen stemming from our previous study are confirmed. The hot X-ray emitting gas 
is shown to have a multi-temperature distribution. 



1. Introduction 

The signature of the star formation (SF) his- 
tory of a galaxy is imprinted in the abundance 
patterns of its stars and gas. Determining the 
abundance of key elements released in the inter- 
stellar medium (ISM) by stars with different mass 
progenitors and hence on different time scales, 
will thus have a strong astrophysical impact in 
drawing the global picture of galaxy formation 
and evolution [8] . It also offers the unique chance 
of directly witnessing the enrichment of the ISM 
[7]. Metals locked into stars give a picture of the 
enrichment just prior to the last burst of SF, while 
the hot gas heated by SNe II explosions and emit- 
ting in the X-rays should trace the enrichment by 
the new generation of stars. 

We have started a project to measure the de- 
tailed chemical abundances in a sample of star- 
burst galaxies, for which we obtained high reso- 
lution infrared (J and H band) spectra with the 
3.6 m Italian Telescopio Nazionale Galileo (TNG) 
and with the ESQ VLT, and both proprietary 
and archival data from the XMM-Newton and 
Chandra missions. Our sample comprises M82, 
NGC253, NGC4449, NGC3256 and the Anten- 
nae, sampling two orders of magnitude in star 
formation, as it ranges from the 0.3 Mq/jt of 
NGC4449 to the - 3OAf0/yr of the Antennae and 
NGC3256. 
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Preliminary results for M82 were achieved with 
the available XMM-Newton archival data, and 
hinted for a confirmation of the expected scenario 
in which the gaseous component has a higher con- 
tent of a-elements than the stellar one, and a sim- 
ilar content of Fe [10]. However, some new issues 
were posed, since we found a very low abundance 
of O and Ne with respect to other a-elements 
(e.g., 0/Mg - 0.2, Ne/Mg - 0.3) in the hot gas 
present in the central (<~ 1 kpc) regions of M82, 
which could not be satisfactorily explained. For 
this reason, we were granted a deeper observation 
of M82. A report on our preliminary analysis ap- 
peared in [12]. Some issues concerning the data 
analysis of complex and well-sampled spectra by 
standard "black box" tools, such as XSPEC, the 
Mekal and APEC models, etc. are presented in 
this contribution. The main paper is currently in 
preparation. 

2. EPIC data reduction at the limits of 
standard tools 

M82 is a bright galaxy with an ongoing star- 
burst and a complex spectrum. Thus, the data 
reduction of such a deep observation (100 ks) 
has proved difhcult and time-consuming, posing 
many challenges about the involved astrophysical 
processes and the available software apparatus. 

The EPIC data were analysed with the SAS 
software, version 6.5.0. After screening for back- 
ground flares, about 73 ks of data were accepted 
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for analysis. We began our analysis from the cen- 
tral region of M82, by extracting spectra from a 
circle with diameter 1' and centred on the coordi- 
nates 09:55:51 and +69:40:39. Background spec- 
tra were extracted from the blank-sky data files 
[13] after normalisation to the background levels 
observed in the M82 data. 

The largest emission from hot plasma is found 
in the centre of M82. However, many point 
sources are also present in the same region. To 
choose a model for the point sources, we anal- 
ysed the spectrum in the 3-8 kcV band, and we 
found that it can be well described by a power 
law with photon index T = l.GOlg Qg. This power 
law has been always included in the models dis- 
cussed below regarding the centre of M82. A Fo 
thermal line is also detected at 6.66 ± 0.02 keV 
with equivalent width ~ 100 eV. 

The abundances for O, Ne, Mg, Si, S, Fe were 
left free to vary, since these elements have strong 
lines in the considered wavelength range. Other 
a-elcmcnts with weaker lines (Na, Al, Ar, Ca), 
which would not be sufficiently constrained on 
their own, were bounded together using solar ra- 
tios, accordingly to the proximity in their atomic 
number, e.g Na and Al to Mg, Ar and Ca. 

2.1. Multi-temperature plasma 

In the literature, the gaseous emission is usu- 
ally modelled with one or two single-temperature 
components. A 'warm' component usually 
aromid 0.7 1 keV is almost ubiquitous [11]. A 
high energy component is also often reported [3,2] 
in studies based on ASCA or BeppoSAX data, 
described cither as a power-law, or as a thermal 
spectrum with kT ~4-10 keV. Because of the 
poor spatial resolution of those satellites it is not 
clear if this component is really diffuse or rather 
due to point sources. Thus, some care should be 
taken in comparing literature results with this pa- 
per, since in the following we'll account for both 
point sources (via the power law discussed in the 
previous paragraph) and for hot gas. A 'cold' 
component, around 0.1 0.3 keV, is also found 
sometimes [4]. In general, the better the quality 
of data, the larger the number of components re- 
quired to provide an acceptable fit. Thus, the use 
of multi-temperature models that allow for a fine 
control of the DEM, and are also insensitive to the 
Fc-bias [1], might be preferable but this is hardly 
feasible, because of demanding requirements in 
data quality and computing time. The grade of 



our M82 observation is sufficiently high to allow 
for the first time the use of multi-temperature 
models even in spatially-resolved spectroscopy. 

The results of the fit showed that a hot com- 
ponent {kT >~ 5 keV) is required in order to 
fit the spectrum. Several versions of the com- 
mon MEKAL code [9] are present with multi- 
temperature flavours in the XSPEC distribution 
(version 11.3.2g), among which the most flexible 
are c6vmekl and c6pvmkl. These models parame- 
terise the plasma DEM, i.e. the amoimt of plasma 
with a given temperature, with a 6"^ order poly- 
nomial, and allow for variable abundances. The 
difference between them is that c6pvmkl restricts 
the DEM to positive values, while c6vmekl also 
allows negative values: this means that negative 
amounts of plasma are permitted. Thus we chose 
to use c6pvmkl. 

However, these models allow only temperatures 
with kT < 7 keV, causing the b(^st-flt DEM to 
have a sharp cut, clearly noticeable in Fig. 1 (left 
panel). Feeling that this limitation was some- 
what artiflcial, we modifled the available multi- 
temperature plasma models by enlarging the allo- 
cated energy range, so that the high-temperature 
part of the DEM could be better sampled. The 
best fit DEM for both the standard and 'enlarged' 
models are shown in Fig. 1. The improvement in 
was however not significant, indicating that 
the DEM sampling was already sufficiently good: 
as it may be seen in Fig. 1, the high- temperature 
component is distributed like a rather narrow bell 
curve, so little difference is found between a sin- 
gle temperature component and the sum of two 
components with very similar temperatures. 

We also checked if the hot component 
could have a non-thermal origin: since star 
forming galaxies are known as radio sources, 
bremsstrahlung emission should in principle arise 
from the same population of electrons from which 
the synchrotron emission is originated. Thus we 
developed a non-thermal bremsstrahlung model 
for XSPEC. We found that the hot thermal com- 
ponent could not be replaced by a non-thermal 
one, because of i) the fit statistics worsens; ii) 
given the observed X-ray flux from the non- 
thermal bremsstrahlung component, and hence 
the required population of free electrons, an un- 
realistically low galactic magnetic field would be 
required to account for the observed radio (syn- 
chrotron) emission. 
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Figure 1. The best fit DEM for the central regions of M82 by using the standard c6pvmkl model (left 
panel) and its enlarged-energy-array version (right panel). The total spectrum is computed by adding 
several components, one for each sampling point shown here. 



2.2. Unidentified lines 

There seems to be two lines in the data which 
are not present in the model (see Fig. 2). If the 
lines are parameterised with Gaussian models, 
the best fit energies are 1.22 ± 0.01 keV (equiva- 
lent width ^ 17 eV) and ~ 0.78 keV (equivalent 
width ~ 20 eV). 

These lines are present through all the spatial 
extent of the outflow, which is about the central 
10' of the MOS and PN chips. By looking at the 
Chianti line emission database [5] the best candi- 
dates for the 1.22 keV line are the Fe XVII line 
at 1.223 keV or the Ne X line at 1.211 keV, while 
for the 0.78 keV line the best candidate is the Fe 
XVIII Hne at 0.7832 keV. Although i) the effec- 
tive areas in all EPIC cameras in the considered 
energy range are a smooth function of the wave- 
length and do not show any feature around the 
considered wavelengths, and ii) no background 
feature is present near the considered wavelengths 
(moreover, the count rate of the central region 
of M82 is about 100 times larger than the back- 
ground), we believe that these lines might also be 
instrumental, background or model features, for: 

• the 1.22 keV line is flanked at longer wave- 
lengths by an absorption-like feature. A 
close look to Fig. 2 suggests that a line, if 

indeed present, should be accompanied by 
problems in modelling the continuum. 

• The 1.22 keV line is not present in the RGS 
data (Fig. 3), and the other one at 0.78 keV 



seems to be present with lesser intensity. 
However, this might also be due to the RGS 
counts being ~ 1/10 of the EPIC ones. 

• The line intensity is reduced if the APEC 
model is used inplace of Mekal (see below) . 

To further check whether the presence of the 
two 'unidentified lines' might be related to the 
details of the spectral model, we fltted the data 
using the APEC model [14]. Since no multi- 
temperature APEC model is present in the cur- 
rent XSPEC distribution, we built one based on 
the Mekal one, with the only difference that the 
APEC routines are called. Switching from the 
Mekal- to the APEC-based model improves the 
fit: for the Mekal one, one gets, considering only 
the MOS spectra, = 1459 with 889 degrees of 
freedom (d.o.f.; x^r = 1.64); for the APEC one, 
= 1274 {x^r = 1.43). A visual inspection of 
the residuals confirms that the APEC model pro- 
vides a better representation of the data (Fig. 3). 
As a technical note, we point out that no signifi- 
cant difference in the value of was obtained for 
the Mekal model with the switch parameter set 
to "calculate" , as opposed to the default value of 
"interpolate from a pre- calculated table" (see the 
XSPEC documentation). 

Thus, we decided to treat these lines as non- 
astrophysical features. A detailed comparison of 
the calculations from the two considered spectral 
codes is being performed, and its results are de- 
ferred to the main paper in preparation. 
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Figure 2. EPIC spectra of the central (2') region of M82. Green points (upper data series): FN. Red and 
black points (lower data series): MOS. Some of the strongest emission lines are indicated. The arrows 
below Fe and O show line blends. The position of the two 'unidentified lines' is also shown with a question 
mark. 



3. RGS data 

The RGS spectra for the inner region were ex- 
tracted with the rgsproc tool, considering only 
events within 90% of the FSF width (i.e. ~ 1') in 
the cross-dispersion direction. The background 
spectra were taken from the blank-sky observa- 
tions using the rgsbkgmodel tool. The line width 
in the RGS data is dominated by the M82 size 
2') so the rgsxsrc XSPEC model, which con- 
volves the instrumental line width with the source 
spatial profile, has been used to account for the 
source extent. This model needed to be patched 
to allow the joint fit of grating and CCD data. 
The spectra are shown in Fig. 3. The best fit 
parameters from fits to either the joint EPIC- 
RGS data and RGS data alone are consistent with 
being a spatially-weighted average of those dis- 
cussed above (see Fig. 5). 

4. Spatially resolved spectroscopy 

We report here on the preliminary analysis con- 
duced on both the northern and the southern out- 
flows, mainly making use of EPIC. In order to 
study the different properties of the hot gas as 
it flows and/or is heated from the central star- 
burst towards the intergalactic space, we divided 
both outflows in five sub-regions. Each region 
has a rectangular shape, the larger side being 



parallel to the galaxy major axis. The spectra 
were extracted from the MOS (0.5-8.0 keV) and 
PN (1.0-9.0 keV) data, and fitted with the same 
model used for the central region, the only differ- 
ence being point sources, which in the outflow are 
sufficiently apart to be excluded when extracting 
spectra. 

The best-fit chemical abundances from EPIC 
and RGS data are shown in Fig. 5 along with 
results from our previous paper [10] relative to 
infrared observations for the central regions. The 
abundances are scaled according to the Grevesse 
& Sauval [6] solar composition. No significant 
changes are found in the temperature of the 
plasma. Only one region required in its spectrum 
some absorption in excess of the Galactic value. 
Our previous finding, that the inner region of M82 
is somewhat devoid of the lighter a-elements, is 
thus confirmed. Moreover, it is found that all ele- 
ments are more concentrated in the outflow than 
in the centre. The centre/outskirt abundance ra- 
tio ranges from ~ 1/10 for O to 1/3 for the 
heavier elements (Mg, Si, S). 

In M82 both the hot gas and the stellar phases 
trace a very similar Fe abundance. Indeed, since 
Fe is mainly produced by type la supernovae 
(SN), it is expected to be significantly released in 
the ISM only after a few 100 Myr from the local 
onset of star formation. At variance, a-elements 
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Figure 3. Residuals of EPIC data with different models. The model parameters are best-fit ones and 
thus are slightly, yet not significantly, different in the two cases. The reduced are 1.64 for Mekal and 
1.43 for APEC. 



(O, Ne, Mg, Si, S) are predominantly released 
by type II SN with massive progenitors on much 
shorter time scales. 

The overall a-elements/Fe enhancement is con- 
sistent with a standard chemical evolution sce- 
nario only for the heavier elements (Mg, Si, S). 
The lighter elements (O, Ne) show a different dis- 
tribution: in the inner regions of the galaxy, the 
0/Fe and Ne/Fe ratios appear to be too low for 
this scenario to hold. In fact, the X-ray derived O 
abundance in the centre is less than the infrared- 
derived one. 

5. Conclusions 

The analysis of a very deep (~ 100 ks) observa- 
tion of the starburst galaxy M82 has opened some 
new questions, which may be summarised as: 

• there is a discrepancy in the oxygen abun- 
dance between the hot gas and stars, which 
goes in the opposite direction with respect 
to the expectation of chemical evolutionary 
models; 

• strong vertical abundance gradients in some 
metals have been discovered; 

• the X-ray emitting gas has a "multi- 
temperature" nature. 

We have worked out some possibilities for in- 
terpreting these results, which will be presented 
in the forthcoming main paper. 
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Figure 4. RGS spectrum of M82. Black, RGSl; red, RGS2. Only the first order spectra are shown for 
simplicity. The residuals add almost all in the O VII and O VIII lines. 




Figure 5. Variation of chemical abundances with increasing height on the galactic plane. Only O and 
Fe+Ni are shown for brevity. Black: abundances from X-ray MOS data. Red: abundances from infrared 
data (corresponding to red supergiant stars in the galaxy central region). Blue: abundances from joint 
EPIC-RGS fits. 



